The attainable superconducting Tc in a model of 
phase coherence by percolation. 
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The onset of macroscopic phase coherence in superconducting cuprates is considered to be deter- 
mined by random percolation between mesoscopic Jahn- Teller pairs, stripes or clusters. The model 
is found to predict the onset of superconductivity near 6% doping, maximum Tc near 15% doping 
and Tc — T* at optimum doping, and accounts for the destruction of superconductivity by Zn dop- 
ing near 7%. The model also predicts a relation between the pairing (pseudogap) energy and Tc in 
terms of experimentally measurable quantities. 
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The occurrence of pairing and a phase coherent super- 
conducting state in the cuprates is thought to occur inde- 
pendently, especially in the so-called underdoped state. 
This is contrary to the behavior in BCS superconductors, 
where they occur simultaneously. This enables us to dis- 
cuss the mechanism for the establishment of phase coher- 
ence in itself, starting with pre- formed real-space pairs. 
In this letter we investigate the possibility of the for- 
mation of a macroscopic superconducting phase-coherent 
state by the spreading of phase coherence by percolation 
(PCP), calculating the attainable Tc and some of the 
critical points in the phase diagram. 

We assume that superconductivity in the cuprates 
arises when pre-formed pairs acquire mutual phase co- 
herence. The pairs are considered to be described by 
the mesoscopic Jahn- Teller (MJT) pairing model, which 
gives a detailed mechanism for pair formation arising 
from degeneracy of hole orbitalsjl, |^. In this model the 
size of each pair Ip (including deformation of the sur- 
rounding lattice and spins) is less than or equal to the 
coherence length but larger than the lattice constant 
a, i.e. CL < Ip £,s- The percolation model proposed here 
is considered primarily for the underdoped state which is 
mixture of Bosonic pairs and Fermionic excitations 
separated by a pairing energy gap Ep ^ . 

Let us now assume that macroscopic phase coher- 
ence is established when the percolation threshold is 
reached and pairs touch each other to form a phase co- 
herent path throughout the crystal. The situation con- 
sidered is (ij/namic|P, where the pair presence fluctu- 
ates statistically on a timescale given approximately by 
T ~ h/ EpCxp{Ep/kT). Figure la) shows a snapshot of 
the Cu02 planes filled by MJT pairs for 6% doping at 
T — (all pairs in the ground state). The average size of 
the pairs is given by the experimentally determined value 
of 7 = 7r/Zp, where 7 is the range of the interaction in k- 
space as given by the inelastic neutron scattering (INS) 
experiments 0. The area occupied by each pair in the 
Cu02 plane is: 



The color (density) depicts the amplitude of the lattice 
deformation caused by the paired carriers predicted by 
the MJT model. The situation depicted in Figure la) 
can be described well by the 2D bond percolation (BP) 
model shown in Figure lb), whereby each possible site 
for a pair is depicted by a dot, and the occupancy of the 
site is depicted by a bond to the next dot. 

The percolation threshold for the BP model in 2D is 
known to occur when the "volume fraction" Fy — 1/2 [8]. 
The "volume fraction" occupied by the pairs at temper- 
ature T and doping density x is then given by: 



Fv = ^xnp{T)Vp/V 



(2) 



(1) 



where np{T) is the thermal pair occupancy at tem- 
perature T and V = a? is the volume of the unit cell, 
and a is the in-plane lattice constant [10]. The factor 1/2 
arises because two doped carriers are required to create 
one pair. 

Let us calculate the critical pair density at the percola- 
tion threshold at T = 0, which corresponds to the onset 
of superconductivity on the underdoped side (point A 
in Fig. 2). In this case all the carriers are paired in the 
ground state and np{T = 0) = 1, so Xc = 2FvV/Vp. 
Assuming Fy — 1/2 and that 7 ~ 7r/2a from the INS 
experiments^, we obtain Xc = 0.08. The experimental 
value for the onset of superconductivity for most cuprates 
is indeed very close, being 0.06 ± 0.01 for La2-2;Srj;Cu04 
|ll| and Xc = 0.05 ± 0.01 for YBCO |l| for example. 
(Note that the insulator-to-metal transition does not nec- 
essarily coincide with this point, since the percolation 
threshold for metallic behaviour may be different p3[ ) . 

Next, let us consider the situation near optimum dop- 
ing. Here, the optimum carrier concentration is typically 
15%. From (2), at T = and x^pt = 0.15 (point B in Fig 
2) we obtain Fy = 0.94. Thus near optimum doping, the 
Cu02 planes are fully occupied by pairs. Further doping 
results in overlapping pairs [ij] and the formation of 
more spatially extended metallic regions 

However, for T > 0, the population of pairs is reduced, 
and at some finite critical temperature Tc, the percola- 
tion threshold is reached (point C in Fig. 2). Using Eq. 2 
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with Fv = 1/2 and the experimental value of 7, we obtain 
a pair occupancy np{T) ~ 0.53. This means that near op- 
timum doping, at Tc the pair occupancy is approximately 
half, which is also the point where kT* ~ Ep (Wc follow 
the convention that the temperature at which the ther- 
mal energy becomes comparable to the pair binding en- 
ergy is called the pseudogap temperature T*). Thus the 
percolation model predicts that Tc ^ T* near optimum 
doping, which is indeed the case in most cupratesp^. 

Considering the population of pairs as a function of 
temperature in the MJT model[|[ we can calculate 
Tc as a function of doping. We start with the expression 
for the number of pairs np{T) at temperature T predicted 
by the MJT model §: 

np{T) ~ 1 - K{Ep/kT)-^/'^ cxp ( - Ep/kT). (3) 

Ep is the pair binding energy, and K in a, constant 
describing the density of unpaired states. Substituting 
into (2) and solving for T, the temperature Tc where 
Fv— 1/2 is approximately given by: 



Tc{x)/Ep^l/HK- 



-)• 



(4) 



To test this relation, in Figure 3a) we plot the predicted 
Tc using experimental values of Ep from quasiparticle re- 
combination data[16] on YBCO as a function of doping 
and using the value oi K — 200 obtained from experi- 
mental fits of the T-dependence of Raman and neutron 
scattering data|^. Remarkably the magnitude of Tc is 
near 90K near optimum doping x = 0.15. A universal 
" pseudogap" ratio of Ep /Tc as a function of doping x 
given by Eq. (4) is plotted in Figure 3b). Over most of 
the doping range it is near 6 and increases near the per- 
colation threshold Here the value of Ep /Tc 
depends only on K and Xc, which can be determined in- 
dependently from many different experiments!^ ^ 
and is in agreement with measured values from single- 
particle tunneling experiments joj. 

In the model above we have not discussed the mecha- 
nism for establishing phase coherence between pairs (or 
clusters of pairs). Since typically the pair size is approx- 
imately equal to the coherence length Ip ~ 7r/27 ~ 
we describe the mechanism in terms of pair Josephson 
tunneling, as was already discussed briefly in ref. [Q. 
The pairs together with their deformation are considered 
immobile on the timescale of h/kTc, which makes per- 
colation relevant . However, the charge carriers can 
tunnel from one pair region to the next in pairs, form- 
ing a phase coherent state throughout the crystal. This 
tunneling occurs resonantly between the Bosonic states 
of each MJT pair-region, without carrying the deforma- 
tion itself. The tunneling amplitude is independent of 
the presence of other particles within the deformed pair- 
region, because the deformed regions are sufficiently large 
to accommodate more than 2 particles. However, if the 
final cluster is populated by a pair, the Coulomb charging 
energy is given by a Hubbard term U '~ 4e^/(eo'p) and 



we may describe the system in terms of a Boson-Hubbard 
model imposed on a pseudo-random bond-percolation 
network. (This is somewhat different from the scenario 
for superconductivity previously considered in the Boson- 
Hubbard model on a regular lattice]!^.) 

The simple percolation model ignores completely the 
effects of thermal phase fluctuations, which are known 
to be important in cuprates and other systems such as 
charge-density wave systems |l^. The effect of fluctua- 
tions is to reduce the Tc, especially in the strongly un- 
derdoped region, where the superfluid density and the 
phase stiffness are small[^. This could explain the dis- 
crepancy between the predicted Tc and the actual Tc for 
small doping (Fig. 3a). 

We can apply the same percolation arguments to the 
discussion of the role of metal impurities such as Ni and 
Zn which replace Cu in the CuO planes. By a modified 
" swiss-cheese" type of argument pl[, if we ignore impu- 
rity aggregation, the real-space area occupied by each 

impurity is of the order of tt . This space is unavail- 
able for the formation of pairs. (The justification for this 
comes from many different experiments, the most direct 
being the scanning tunneling microscope measurements 
of Pan et al (2000)|Q and muon spin rclaxation|21 .) By 
equation (2) the critical concentration of impurities for 
the total destruction of superconductivity (i.e. percola- 
tion threshold at T = 0) is xf " = 8%, close to the 7% 
typically found experimentally for Znpl|. Of course dif- 
ferent ions may inhibit the formation of mesoscopic pairs 
over slightly different distances and the exact value of the 
critical impurity concentration will vary. 

The effect of increasing carrier density is to make ai 
gregation of pairs into stripes more and more probable 

Even at the percolation threshold near 6% doping 
Fig. la), the planes are substantially occupied and it is 
natural that the formation of longer stripes will occur, 
especially at low temperatures. The additional energy 
gain associated with the formation of such stripes will be 
mainly associated with the sharing of the deformation 
and the energy scale for the aggregation is expected to 
be significantly smaller than the pairing energy (pseudo- 
gap). Consequently the formation of stripes (such as is 
shown in Fig. la) will occur at temperatures below T*, 
and their density will be determined by statistics. How- 
ever, if these larger stripes or clusters are metallic, they 
occupy space, but do not contribute to the superconduct- 
ing volume fraction. Although superconductivity can be 
induced in these regions by the Bosonic pairs via the 
proximity effect, increasing doping will lead to a reduc- 
tion in the volume occupied by pairs and hence a rapid 
decrease in Tc- The drop of Tc in the overdoped phase 
may then be described by the same PCP mechanism as- 
suming the appearance of metallic stripes results in a 
reduced volume fraction of pairs and eventual disappear- 
ance of superconductivity at point D in Fig. 2. 

In spite of the relative simplicity of the PCP model, 
with very few underlying assumptions we are able to re- 
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produce many features of the cuprate phase diagram ir- 
respective of the symmetry of the pairing channel. The 
model gives a theoretical prediction for the "pseudogap" 
ratio Ep/Tc in terms of experimentally measurable quan- 
tities. It enables us to quantitatively predict the critical 
carrier concentration for the occurrence of superconduc- 
tivity on the underdoped side as well as estimate the 
from experimentally measured values of the pairing en- 
ergy (or pseudogap) Ep. 

One of us (VVK) acknoweledge the financial support 
from EPSRC (UK, Grant R46977). 

Figure 1. a) the amplitude of the lattice deformation 



caused by pairs described by the mesoscopic Jahn- Teller 
model. The picture corresponds to a "snapshot" at 6% 
doping at T = 0. b) the bond-percolation model describ- 
ing the situation in a). 

Figure 2. A schematic phase diagram of the cuprates 
with points referred to in the text. The dashed line rep- 
resents the pseudogap temperature T* . 

Figure 3. a) as a function of x obtained from mea- 
sured values of Ep in YBCO. The dashed line approx- 
imately represents the measured T^. b) A plot of the 
pseudogap ratio Ep/kTc predicted by eeq. (4). 
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